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Abstract. For a large class of lattice models we study the ther- 
modynamic factor, $, of the collective surface diffusion coefficient 
near a first-order phase transition between two phases at low tem- 
peratures. In a two-phase regime its dependence on the coverage, 
9, is <E> « 0/[{0 - 0-)(0+ — 9)N], where N is the number of adsorp- 
tion sites and 6± are the single-phase coverages at the transition. 
In the crossover between the two-phase and single-phase regimes 
&(0) is shown to have a more complex behavior. The results are 
applied to a simple 2D lattice model. 



Surface diffusion is an intriguing phenomenon that has attracted 
much attention of surface scientists in the past few decades [H12]- One 
of the relevant transport coefficients for surface diffusion is the collec- 
tive (or chemical) surface diffusion coefficient, D c . It is associated with 
the decay of fluctuations in the adparticle density at large time and 
space intervals and is defined via the Fick's first law, J = —D C V9, 
where J is the surface diffusion flux and 9 is the surface coverage. 

Considering suitable 2D lattice gas models, computer simulation 
methods have been frequently applied in theoretical studies of D c (and 
of surface diffusion in general) with a particular interest in the presence 
of phase transitions and their effects on surface diffusion [3], HI El El [7J 
M, M ELD H2]- Usually, the observation that D c = $ D CM [lj was used, 
where -Dcm is the center of mass diffusion coefficient and represents ki- 
netic properties of D c , while $ = (N oc )/((N% c ) — (N oc ) 2 ) is associated 
with thermodynamic properties of D c and is called the thermodynamic 
factor (iV oc is the number of adsorption sites in a system occupied by 
adparticles and (•) denotes the statistical mean value). The factor is 
simply related to the coverage 9 = (N oc )/N and isothermal suscepti- 
bihty X = /3((iV 2 c )-(iV 0C )2)/iV ) 

$ = (1) 

X 
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where N is the total number of adsorption sites in the system and 
ft = l/ksT is the inverse temperature. In the hydrodynamic limit 
(when 9 varies slowly with space and time) D c can be obtained from 
the system free energy, and so it is a purely thermodynamic quantity 

nana lis]. 

In this letter we use a rigorous statistical mechanical approach to 
study $ near a first-order phase transition with a (in general asym- 
metric) coexistence of two phases (denoted as '+' and '— '). To this 
end, we consider a finite array of N adsorption sites with periodic 
boundary conditions and use that, for a wide range of classical lattice 
spin models (such as the models with a finite range m-potential and 
a finite number of ground states), the coverage and susceptibility in 
the array have a remarkably universal behavior near such a transition 
[TBI [T7] . Hence, the behavior of $ near the transition as a function of 
the chemical potential, /i, as well as of the coverage 6 can be described 
in a very general and unifying way. It is only needed that the temper- 
ature is low (P is large) and that the array is large and not too oblong 
(N is large and the number of boundary sites in the array is of order 
of a fractional power of N). 

The [i- dependence o/$ near the transition. Generally, the coverage 
can be expressed as [16] 

9 = (-/;)A+ + (-/i)A_ + R N , (2) 

where f± and — f± are single-phase specific free energies and coverages, 
respectively, the weight factors A± = 1/{1 + exp[±/3(/ + — f_)N]}, and 
the error term R N = 0(exp(— cfiyN)) with c > (primes denote 
derivatives with respect to /i, and 0(x) stands for a term that can 
be bounded by const a;). To get an explicit /^-dependence of 6 near a 
transition point, // t , one applies Taylor expansions around fi t , yielding 

d = Y + — tanh 2 + 0(/?2|/i ~ ^ |2) + Rn (3) 

(see Fig. EEJa)). The shorthands qo = q+ + q~ and Aq = q + — g_ denote 
the sum and difference, respectively, of any single-phase quantities q±, 
#± = 0± + X±(j*-fh), and Y = f3N{(i- ii t )[A6 + A X (//-// t )/2], where 
0± = —f±(jh) an d X± — ~f±(f't) are the single-phase coverages and 
susceptibilities at the transition. If the plus (minus) corresponds to the 
phase stable above (below) fi t , the coverage jump at \i t is Q(p, t + 0) — 
9(fi t — 0) = AO > 0. For \i beyond the transition Eq. (T5]) yields another 
expression for 9 [16] . 

9 = 6 00 + P N + R N . (4) 
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The infinite- volume coverage 6*00 is equal to —f' + (—fL) for \i above 
(below) fx t and P N = 0(exp(— C\fi — fM t \(3N)) with C > 0. Using 
s = sign(/i — fi t ) both for numbers ±1 and for subscripts '±', we may 
write Boa = -ft. 

Equation (j2J) holds also for the /i-derivatives of 9 with the same error 
term [16]. So, we can readily obtain a general expression for the 
susceptibility, 

X = (-Af) 2 (3N\ + \_ + (-f';)X + + (-/^)A_ + R N , (5) 

where — f± are the single-phase susceptibilities. The /i-dependence of 
X then follows analogously to that of 9: either from Taylor expansions 
around fi t (useful near the transition), 

Al 2 /3iV xo Ay Yl 

x = w?x + T + — tanh y (6) 

+ 0(/3 3 iV|/i - /i t | 2 ) + 0(/3 2 |/i - fi t \) + R N (7) 

with Yi = A9(3N(fi—fi t ) (see Fig.[T](b)), or by using the infinite-volume 
susceptibility Xoo = —f'J (useful beyond the transition), 

X = Xoo + 0((3NP N ) + R N . (8) 

These formulas for 9 and x immediately yield the /i-dependence of 
the thermodynamic factor $. Near the transition region, — fi t \ <C 1, 
it is given by Eqs. ([3]) and ([6]) (see Fig. [H(c)), while farther away from 
the transition, — [i t \ ^> 1/N, it is given by Eqs. PJ and (jSJ). 

At a certain distance from \i t (both below and above fit) there are two 
narrow regions of \x where 9 and x reduce to their infinite- volume limits 
and Xoo so that the behavior of $ switches from a two-phase regime 
to one of the two single-phase regimes (see Fig. [1(c)). According to 
Eqs. ([3]) and (Q, these two crossover regions occur when {3N exp(—\Y\) 
becomes small, i.e., within the intervals 5fi- ^ |/i — fit\ ^ 5fi+, where 
5fj,± = (InN ± K\n\nN)/A9f3N with K > large. 

The 9 -dependence 0/$ near the transition. The coverage dependence 
of $ follows upon the elimination of \x between 9(\x) and This 
procedure is simple to achieve within the two-phase regime, |/i — /i t \ ^ 
8fi-, because it reduces to the elimination of Y between Eqs. fl3]) and 
dSD- We readily get 

Q 

$= (9-9_)(9 + -9)N + 0(\nN) ' ^ 



where t± = #(/if ± 8^-) are the most extreme values that 9 attains in 
the two-phase regime. Since t± = 0± T A9\n K N/N + O (In N/N) by 
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Figure 1 . The dependencies of (a) coverage, (b) suscep- 
tibility, and (c) thermodynamic factor on the chemical 
potential near the transition (full lines). For two larger 
sizes of a system the dependencies are depicted by dashed 
and dotted curves. The points Xmax ~ (A9/2) 2 /3N and 
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Figure 2. The dependence of $ (the full line) in the 
interval £_ ^ 9 ^ t + (the two-phase regime) as given by 
Eq. OH]). For two larger sizes of a system the dependence 
is depicted by a dashed and a dotted curve. A mini- 
mum of $ is attained at 9 m - m ~ (6_6 + ) 1 / 2 and its value 

iS $ min ^ l/(y/r-y/9^) 2 N. 

Eq. 03]), we have 6L < £_ < t + < 9 + . The behavior of $ as given by 
Eq. (Q is shown in Fig. [2] 

Next, consider the range <5/i_ ^ — fi t \ ^ dfx with d/i = 1/ A9f3N 3 / A 
that contains both crossover regions as well as a part of either single- 
phase regime (instead of N 3 ^ any power N a with 1/2 < a < 1 can be 
considered). To eliminate \i between 6*(/i) and x(jt/) we rewrite Eqs. fl3]) 
and ^ as 6» = S + - sAM + Qjy and % = A9 2 (3NA + x s + 

(3NQ N , respectively, with A = exp(-|Yi|) and Q N = 0(hi K N /N 3 ' 2 ). 
Eliminating A between these equations, we get 

A9/3N fn 



with 



uj± = W\ exp 

^ X± 



X± 

where W(y) is the Lambert PF-function (i.e., the inverse toy — W exp W) 
and the restrictions on 9 in Eq. ( flOi) correspond to the two intervals 
of values attained by 9 in the considered //-region below and above fit- 
Thus, r ± = 8(ih ± dfi) = 9±± x±/A9(3N 3 / 4 + 0{l/N 3 / 2 ) by Eq. © so 
that r_ < 9- < t- < t+ < 9 + < r+. 

The behavior of $ as given by Eq. ([10]) is shown in Fig. [3] If 
w± > 1, the two-phase behavior of $ prevails, while a single-phase 
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(a) 




Coverage 



Figure 3. The ^-dependence of <3> (full lines) in the 
crossover regions (a) near Q_ and (b) near 9 + as given 
by Eq. (fTOj) (the error term 0(NQ N ) is neglected). For 
two larger sizes of a system the dependencies are depicted 
by dashed and dotted curves. 

behavior is dominant for u± -C 1. The crossover occurs when u± 
is comparable to 1, i.e., when 9 is comparable to the value, 9* ± = 
9± ± (x±/&6{3N) [ln(A9 2 (3N/x±) - 1], at which u± = 1. For 9 < 9*_ 
(for 9 > 9* + ) the thermodynamic factor is increasing in approximately 
a linear way as (39/x- (as 09/x+), while as 9 crosses 9*_ {9* + ) and gets 
above 6L (below 9 + ) it approaches the two-phase behavior described 
by Eq. ©. 

An example. To illustrate the above results, we consider the following 
one-component lattice-gas model of diffusion on a triangular lattice 
[3j |5l [6]. Each lattice site is either occupied by an adparticle or is 
vacant. The model Hamiltonian is H = sN^ c — fiN oc , where N% c is the 
number of nearest-neighbor bonds occupied by adparticles and e < is 
an attractive interaction between two closest adparticles. The model is 



COLLECTIVE SURFACE DIFFUSION NEAR A PHASE TRANSITION 7 




Coverage 



Figure 4. The ^-dependence of $ for the example 
model with iV = 50 x 50 and 0\e\ = 1.2 In 3 in (a) the 
two-phase regime, (b) the crossover region below fi t , and 
(c) the crossover region above /x t . Circles are values of $ 
obtained numerically upon elimination of \x between 6(fi) 
and x(^) using Eqs. ([2]) and §5§. Full lines correspond 
to formulas (case (a)) and (TTUl) (cases (b) and (c)). 
Dotted lines in (b) and (c) correspond to the two-phase 
regime formula ()9]). 
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equivalent to the standard Ising model, and at low temperatures (for 
/3\e\ > In 3) it exhibits a first-order phase transition at fi t = 3e between 
a fully vacant phase (stable below fi t and thus corresponding to — ) and 
a fully occupied phase (stable above fi t and corresponding to +) [18J. 

As follows from Ref. [16], one has f± = e±+s±, where e+ = 3e— fi and 
e_ = are the specific ground state energies and s± represent thermal 
excitations over the ground states and can be expressed by cluster ex- 
pansions. In the lowest order approximation s± ~ — exp(— f3AH±)/f3, 
where AH + = fi — 6e and API = — /x are the energy excesses of the 
lowest excited states (one adparticle taken from and added to the cor- 
responding ground state). So, —f' + ~ 1 — r]exp[—f3(fi — fi t )], —f'_ ps 
rjexp[/3(fi - fi t )], and -f± « /3r?exp[=F/3(/i - fit)] with r] = exp(/3/i t ). 
The thermodynamic factor <£> of the model can be readily obtained 
from these expressions and Eqs. (j2J) and (jSJ), and it is shown in Fig. HI 
Moreover, in the single-phase regimes, \fi — fi t \ > dfi, Eqs. (jl]) and 
OH]) yield $ w 1 + 0(NP N ) for fi below fi t (i.e., for 6 near 0) and 
$ w ^/(1-0)+O(A^Pjv) for /i above ^ (i.e., 5 near 1). Thus, D c w D CM 
at low coverages, while the Langmuir relation holds at high coverages. 

Finally, D c = D®(3 exp( / 9/i)P 00 /x in the hydrodynamic limit [S], 
where D° is the diffusion coefficient of non-interacting adparticles and 
Poo = 1 — 26* + (df /de)/3 is the probability that two nearest-neighbor 
sites are empty (/ = —ksTXnZ '/N is the system specific free en- 
ergy). Using the results of Ref. [16], we get df/de = [(df + /de)X + + 
(df-/de)X-](l + Rn)- Along with the above formulas, we find that 
at fi = fi t , say, one has ln(P> c /P>°) w -3\e\/k B T - ln[(l - 2rf) 2 N/2] + 
0(1/N). 
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